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Ising lines: natural topological defects within chiral ferroelectric domain walls 

V. Stepkova®^, P. Marton®^, and J. Hlinka^EI 

^Institute of Physics, The Czech Academy of Sciences, Na Slovance 2, 18221 Praha 8, Czech Republic 

(Dated: April 28, 2015) 

Phase-field simulations demonstrate that the polarization order-parameter field in the Ginzburg- 
Landau-Devonshire model of rhombohedral ferroelectric BaTiOs allows for an interesting linear 
defect, stable under simple periodic boundary conditions. This linear defect, termed here as Ising 
line, can be described as about 2 nm thick intrinsic paraelectric nanorod acting as a highly mobile 
borderline between finite portions of Bloch-like domain walls of the opposite helicity. These Ising 
lines play the role of domain boundaries associated with the Ising-to-Bloch domain wall phase 
transition. 

PACS numbers: 77.80.-e,77.80.Dj,77.84.-s 


Perovskite ferroelectrics are the key materials in the 
current ferroelectric research and a considerable effort 
was recently devoted to the investigations of their do¬ 
main structure and their domain wall properties. Much 
less is known about ferroelectric line defects, although 
this topic is also gaining attention [l|-l5|, It hIII l33| . for 
example, in relation to magnetic vortices and skyrmions 
UM- Ferroelecric domains and line defects may be in¬ 
spiring in other research areas, too, for example in the 
particle physics |12 |. 

One of the most intriguing recent result in this field 
is the prediction of so-called Bloch-like domain walls 
in BaTiOs and PbTiOs, based on Ginzburg-Landau- 
Devonshire model [l3|-[l6| and first-principles calcula¬ 
tions [i30- Previously, we and others studied [211]- 
oriented 180-degree domain walls in the rhombohedral 
BaTiOs, separating domains with spontaneous polariza¬ 
tion parallel and antiparallel to [111] (as usual we re¬ 
fer to the cubic axes of the parent paraelectric phase of 
BaTiOs, see Fig.[T]a). Calculations carried out within 
the Ginzburg-Landau-Devonshire model of Refs.|l9|, I2 qI 
suggest that this domain wall has a Bloch-like (bistable 
and chiral) structure 1110. At the same time, it has 
been shown that this ([211]-oriented) Bloch-like wall can 
be transformed to an achiral, Ising-like domain wall by 
a moderate uniaxial stress [I 3 • The transformation pro¬ 
ceeds as a continuous, symmetry-breaking phase tran¬ 
sition associated with the divergence of dielectric per¬ 
mittivity [2l| and with disappearance of the polarization 
Pdw within the domain wall interior [l5[. The internal 
domain wall polarization Pdw, existing only up to the 
Bloch-to-Ising transition point, is parallel or antiparal¬ 
lel to [oil] direction. It is natural to ask whether the 
antiparallel Pdw states could coexist within the same 
domain wall plane, similarly as the ferroelectric domains 
may coexist in the bulk of ferroelectric crystal 


In our previous work 0 we have studied the case 
when the Pdw vector has an opposite orientation within 
neighboring domain walls, as indicated in Fig.[T]b. Such 


domain structure contains only the domain walls of the 
same helicity and it is thus a chiral structure as a 
whole. In contrast, in the present work we have ex¬ 
plored a racemic domain structure, which contains an 
equal amount of left-handed and right-handed areas on 
average as well as within each domain wall plane itself. 
This arrangement, depicted in the Fig.[T]c comprises an 
interesting topological defect - a borderline between do¬ 
main wall regions of opposite helicity. This defect could 
be in general denoted as a polarization disclination line, 
recalling for example the analogy to the dislocation lines 
of crystal lattices and disclination lines of liquid crystal 
textures [25[. However, for the reasons given below, we 
prefer to denote the specific defect studied here as an 
Ising line. The aim of the present paper is to describe 
the basic properties of this so far unexplored structural 
object. 

In order to learn about the structure of the Ising 
line, we have carried out standard phase-field sim¬ 
ulations based on numerical solution of the time- 
dependent Ginzburg-Landau equation for a phenomeno¬ 
logical Ginzburg-Landau-Devonshire energy functional 
complemented by exact calculation of the long-range elec¬ 
trostatic energy associated with the inhomogeneous pro¬ 
file of the polarization field [ 2 ^. The local mechani¬ 
cal equilibrium and strain compatibility is implicitly im¬ 
posed by the method [ 2 ^. Further details about the 
phase-field approach in ferroelectrics can be found in 
Refs.0, [soHai 

In the results displayed here, we employed the same 
Ginzburg-Landau-Devonshire model for BaTiOs as in 
Ref.[l5|. The adjustable temperature parameter of the 
model was set to 118 K, but similar results were ob¬ 
tained at other temperatures within the stability range 
of the rhombohedral ferroelectric phase as well. Peri¬ 
odic boundary conditions were applied to both polariza¬ 
tion and strain fields. Furthermore, the overall stress-free 
conditions were assumed (stress averaged over the peri¬ 
odic supercell is vanishing). The calculation was typically 
conducted at a 128 x 128 x 128 sized discrete mesh with 
0.5 nm or 0.25 nm spatial steps. This allows to inspect 
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FIG. 1: Predicted domain structure with Ising lines. Panel 
(a) shows schematically a domain structure with (211) ori¬ 
ented domain walls between [111] and [Ill] oriented ferroelec¬ 
tric domains. In case of the rhombohedral phase of BaTiOs, 
these domain walls become polar with the inner polarization 
pointing around [Oil] or [Oil] directions. When the inner po¬ 
larization of the wall is alternating from the wall to wall as in 
panel (b), the structure is chiral. Panel (c) shows a racemic 
domain structure with the investigated Ising lines, forming 
the borderlines between areas of the opposite inner polariza¬ 
tion. Orientation of the inner polarization within selected 
domain wall is sketched in panel (d). Panels (e) and (f) show 
the symmetry adapted polarization coordinates and the core 
structure of the Ising lines, respectively. 
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FIG. 2: Profiles of Pr,Ft, and Ps polarization components 
along the [100] direction as obtained from phase-field simu¬ 
lations for rhombohedral BaTiOs at 118 K using the model 
described in the text. Data are taken from the front face of 
the simulation box (a) at the trajectory ”1” shown in Fig.[T]:, 
(b) at the trajectory ”2” shown in Fig.[Tl:, and (c) at the in¬ 
termediate trajectory passing through the Ising lines. 


local polarization vectors at the scale of the perovskite 
elementary unit cell parameter. 

The initial configuration was close to the final relaxed 
configuration shown in Fig.[T]c. A careful design of the 
initial configuration was an essential condition for keep¬ 
ing the Ising lines in the relaxed domain structure, be¬ 
cause these lines are highly mobile and easily annihilate. 
The architecture of our initial configuration obviously re¬ 
spected the basic requirements such as that the energeti¬ 
cally very costly head-to-head or tail-to-tail polarization 
configurations should be avoided. Among others, this 
implies that the borderlines between domain wall regions 
with opposite Pdw are all parallel to Pdw, see Fig.[T]d. 

Subsequently, the initial polarization field configu¬ 
ration has been fully relaxed accordiM to the time- 
dependent Ginzburg-Landau equation [l^, l20|, l26| . The 


final relaxed domain structure state indeed kept the over¬ 
all architecture of the engineered initial state (Fig.[T]c), 
at least as long as the distances between the Ising lines 
were larger than about 3nm and the distance between 
the domain walls was larger than about 6 nm. 


The local arrangement of the local polarization around 
the core of the Ising lines is shown in Fig.[T]f. For a 
more quantitative analysis of the resulting polarization 
profiles, it is useful to employ the symmetry-adapted 
coordinates associated with the orthogonal unit vectors 
r II [111] s II [211], and t || [Oil], see Fig.[T]c,e. For exam¬ 
ple, the Bloch-like nature of the domain walls in the re¬ 
laxed structure becomes apparent when the Pr, Pt, Ps po- 
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FIG. 3: Quasistatic hysteresis loop showing partial switch¬ 
ing of the spontaneous domain wall polarization Pdw by 
the electric held applied along t || [Oil] direction. Results 
were obtained from phase-held simulations at 118 K using 
128x128x128 simulation box with a spatial step 0.5 nm. The 
polarization prohles taken at the points marked A and B are 
shown in Fig. 4. 


larization components are traced across the domain wall 
(see Fig.[2]a). In the heart of the domain wall, where the 
bulk spontaneous polarization changes between [111] and 
[III] orientations, the principal polarization component 
Pr passes through zero. However, at the same time, the 
additional polarization component Pt appears, so that 
Pdw, tbe total polarization at Pr = 0 point, become 
nonzero. The exchange between Pr and Pt components 
on the path across the domain wall corresponds to the ro¬ 
tation of the polarization vector within the plane parallel 
to the domain wall. This is similar to the case of ferro¬ 
magnetic Bloch domain walls 111 113. In the domain 
walls shown in Fig.[2]b, the polarization rotates in the 
same sense as the windings on a right-handed screw go; 
the opposite sense of rotation is encountered in domain 
walls depicted in Fig.[2]a. 

Fig.[2]c displays the polarization profile along a trajec¬ 
tory that connects the antiparallel bulk domain states 
through the Ising line. It demonstrates that within the 
Ising line, the polarization vector vanishes. In other 
words, this profile actually reminds of an Ising wall, and 
this is why the line observed here is denoted as Ising 
line. This strongly contrasts with analogical line defects 
known from ferromagnetism, so-called Neel lines^ where 
the magnetization does not vanish in the domain wall 
center (instead, it only tilts out of the plane of the domain 
wall) [251 . In this sense, the Ising lines obtained here are 
the lines of true singularities of the polarization distribu¬ 
tion and in principle, the direction of the polarization is 
not defined in the core of the line at all. The reason for 
this difference is obviously the importance of the electro¬ 
statics; ferroelectric Neel line would be strongly penalized 
by the associated local depolarization field [39|. However, 
this penalty is less prohibitive than that of the 2D Neel 
wall, and in principle, the ferroelectric Neel line may also 
exist in some ferroelectric materials. 


We have also explored the stability of the domain 
structure with respect to applied bias electric fields. The 
width of the Ising line is about 2nm and such a broad 
thickness ensures its high mobility within the domain 
wall. Concerted drift of Ising lines can be forced by a 
homogenous bias field applied along t || [Oil] direction. 
For example, in the simulations with 0.5 nm step the co¬ 
ercive field sufficient to remove all the Ising lines at 118 K 
was only about 3 kV/mm. 

A hysteresis loop showing the variations of the Pt po¬ 
larization component averaged over the whole supercell 
under a periodically oscillating electric field is shown in 
Fig.[3l This hysteresis is related to the growth of the 
area with a favorably oriented Pt component, mediated 
by the motion of the Ising lines. The comparison of the 
virgin and partially switched polarization profiles within 
the domain wall plane is shown in Fig. 01 The profile of 
Pt component clarifies that the Ising line can move as a 
soliton, without any essential change in its internal form. 

Fig. 0] also shows slight spacial oscillations of the pro¬ 
files of Pg and Pr components. They are partially related 
to a slight bending of the domain wall plane around the 
Ising line (the wall is not completely planar). On the 
other hand, the polarization profiles of the domain wall in 
the fiat regions distant from the Ising line are very similar 
to the ideal polarization profiles obtained in Refs.[l5|, [l6|. 
In spite of a strong electrostatic interaction term in the 
model, domain wall profiles in these fiat regions (Fig.[2lib) 
show a clear indication of the characteristic double-kink 
in the normal (Neel) component of the polarization (Ps). 
It is worth noting that in some high-symmetry domains 
walls, such double kinks on Pg or Pt appear only 

after considering fiexoelectric coupling terms, while 
here, in the case of the [211]-oriented Bloch (bistable [38|) 
wall, the fiexoelectric coupling has actuall y on ly a negli¬ 
gible influence on the domain wall profile [l^, [l6| . In ei¬ 
ther case, the amplitude of these double-kinks (about one 
percent of the spontaneous polarization) is much smaller 
than the magnitude of the Pg oscillations of Fig. 01 intro¬ 
duced by the presence of Ising lines themselves. 


As already mentioned, the Ising lines are parallel to 
the t II [Oil] direction, and in fact, in our simulation the 
polarization happens to be strictly constant along t di¬ 
rection. Therefore, it is convenient to inspect section of 
the simulated domain structure by a plane perpendicular 
to t II [Oil] direction, highlighted in Fig. 5a. Polarization 
distribution in this plane (Fig. 5b) shows that the Ising 
lines in the domain structure investigated here are actu¬ 
ally forming a slightly distorted triangular lattice with 
mutual distances of about 14 nm. 

If we follow the polarization round an oriented circuit 
far from the core of given Ising line (see Fig. 5c), it is 
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FIG. 4: Profiles of Pr, Pt and Ps polarization components 
along r || [111] direction, passing through the Ising line sin¬ 
gularities. Panel (a) and (b) correspond to the points marked 
A and B in Fig. 3 and represent the virgin state and partially 
switched state, respectively. 

found that the polarization rotates once around the s 
direction [28|. Therefore, the Ising line can be consid¬ 
ered as an edge disclination of a unit strength [25|. The 
nearest neighbor Ising lines within the same domain wall 
have always the opposite sense of the polarization rota¬ 
tion. It is worth to note that in case of the flux closure 
vortexes observed in ferroelectric thin films and ferro¬ 
electric nanodots, a similar construction would result in 
rotation around the core of the vortex and would be thus 
analogical to a serew disclination of a unit strength. 

The macroscopic (average) symmetry of the investi¬ 
gated domain structure is mmm, with principal axes 
along the r, s, and t directions. However, all the true 
mirror planes have fractional translation parts. There is 
no simple mirror symmetry plane passing through the in¬ 
dividual Ising line, which is clearly a chiral object. The 
Ising line is translationally invariant along the t direc¬ 
tion. In the continuum limit, the core of the Ising line 
is a line of zero polarization. The Jacobian matrix eval¬ 
uated there has two complex conjugate eigenvalues with 
a positive real part and one zero eigenvalue. This im¬ 
plies that it is a line of singular points of saddle-node 
bifurcation type [^ . 


In summary, we clarified the stability and nature of 
an interesting type of line defect expected to exist in fer¬ 
roelectric Bloch walls as a borderline between areas of 





FIG. 5: Galculated polarization field within the plane per¬ 
pendicular to the Ising lines (a) and its smaller part with 
indicated triangular lattice of the Ising lines (b). Panel (c) 
displays vicinity of one Ising line highlighted in (b) and an 
oriented closed circuit around its core. 


opposite helicity. It is found that this defect is highly mo¬ 
bile and its motion mediates the switching of the Bloch 
domain wall polarization and its helicity. It was found 
that in contrast to Bloch lines and Neel lines, which play 
a similar role in ferromagnets, the core of the ferroelec¬ 
tric defect revealed here is unpolarized. In this sense, it 
is appropriate to term this defect as Ising line. 

We have proposed a simple ferroelectric domain struc¬ 
ture with a triangular lattice of such Ising lines. This 
structure was demonstrated to be stable within the es¬ 
tablished Ginzburg-Landau-Devonshire model for rhom- 
bohedral phase of BaTiOs. It was found that the internal 
structure of the Ising line can be described as a line of 
saddle-node type singular points, forming a polarization 
disclination of a unit strength. In a more simple words, 
this Ising line is a straight and about 2 nm thick para- 
electric nanorod of intrinsic nature, embedded within a 
polar but uncharged 180-degree Bloch ferroelectric wall 
and extended along a {011} type direction. 

The polarization vector within a circuit drown around 
the axis of the Ising line studies here makes a full turn 
around the Bloch wall normal s, so that the defect struc¬ 
ture differs from more frequent cases where the polar¬ 
ization rotates around the axis of the defect, as for ex¬ 
ample in 2D skyrmions or the ferroelectric flux-closure 
vortexes. Similar Ising lines are obviously expected to 
exists in other kinds of ferroelectric Bloch walls. 

This work is supported by the Czech Science Founda¬ 
tion (Project CSF 15-04121S). 
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